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Abstract. A spherical model of skeleton with junctions is investigated 
by Monte Carlo simulations. The model is governed by one-dimensional 
bending energy. The results indicate that the model undergoes a first- 
order transition separating the smooth phase from the crumpled phase. 
The existence of phase transition indicates that junctions play a non- 
trivial role in the transition. 

1 Introduction 

Surface model of Helfrich, Polyakov and Kleinert (HPK) |ll2l,'i| has long been 
investigated for a model of biological membranes [4I5I6I7I813] . A curvature energy 
called the bending energy is assumed in the Hamiltonian, which is discretized on 
triangulated surfaces in numerical studies |1U| . It was reported that the model 
undergoes a first-order transition between the smooth phase and a crumpled 
(or wrinkled) phase |llll2ll3| . Experimentally, a similar transition was recently 
observed in an artificial membrane |14) . 

The curvature energy in HPK model is a two-dimensional one, and it reflects 
a homogeneous structure of membranes as a two-dimensional surface. In fact, the 
conventional picture of biological membranes is connected to such homogeneity, 
and moreover many artificial membranes are known as homogeneous. 

However, membrane skeleton, called the cytoskeleton in biological mem- 
branes, has been considered to play a crucial role in maintaining shape and 
motion. Recent experimental studies reveal that free diffusion of lipids suffers 
from the existence of cytoskeletons: the so-called hop-diffusion was actually ob- 
served |15| . The artificial membrane, in the above-mentioned experimental study, 
is considered not to be a homogeneous surface, because it is partly polymerized 

Therefore, it is interesting to study whether the phase transition occurs in a 
skeleton model, which is obviously different from the conventional surface model 
such as HPK model. The skeleton model is significantly simplified as a model for 
surfaces; it contains only skeletons and the junctions. Consequently, we consider 
that the simplified model enable us to focus on whether skeletons play a crucial 
role on the phase transition in biological and artificial membranes. 

In this paper, we report numerical evidence that a phase transition occurs in a 
skeleton model, which is governed by one-dimensional bending energy. This result 
is remarkable, because no phase transition can be observed in one-dimensional 
object described by Hamiltonian of local interactions. 
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2 Model and Monte Carlo technique 

By dividing every edge of the icosahedron into some pieces of the uniform length, 
we have a triangulated surface. The compartment structure is built on the surface 
by eliminating the vertices inside the compartment. The boundary of the com- 
partment forms the skeletons with junctions. We denote by (N, N$, Nj, L) the 
total number of vertices (including those in the junctions and their neighboring 
vertices), the total number of vertices in the skeletons (excluding both vertices 
in the junctions and their neighboring vertices), the total number of junctions, 
and the number of bonds between two neighboring junctions in a compartment, 
respectively. Note that L — 1 is the number of vertices in the skeleton between 
the junctions. 

A starting configuration of Monte Carlo (MC) simulation is shown in Fig. 
QJa), db), where (N,N s ,Nj,L) = (2562,1440,162,4) and (N, Ns, Nj, L) = 
(19362,10890,1212,4). The dots are clearly seen in (a) and unclear in (b) be- 
cause the size of (b) is relatively large. The surface of the sphere is shown in the 
snapshot for clarity. The hexagon and the pentagon (which are clearly seen in 
(a)) correspond to the junctions. Total number of the pentagon is 12, and the 
remaining junctions are hexagonal. 




(a) (2562, 1440, 162, 4) (b) (19362, 10890, 1212, 4) 



Fig. 1. Starting configuration of surfaces of (a) (N,Ns,Nj,L~) = 
(2562,1440,162,4) and (b) (N,N s ,Nj,L) = (19362,10890,1212,4). Thick 
lines denote the skeletons and the junctions, and small dots the vertices. The 
dots are clearly seen in (a) and unclear in (b) because the surface size is 
relatively large. The surface of the sphere is shown for clarity in both (a) and 
(b). 

Henceforth, we use terminology junction for the hexagonal (or pentagonal) 
objects or for the central vertices of these objects, skeleton for linear object 
between junctions, and surface for junction + skeleton. 
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The Gaussian bond potential Si, the one-dimensional bending energy S 2 

(2) 

for skeletons, and the two-dimensional bending energy S 2 ' for junctions, are 
given by 

si = £(*i-*» 2 . s { 2 l} = j2a-u-t J ), sf = £(1-^.11,.), a) 

to) to] to) 

where m i s the sum over bonds (ij) connecting all the vertices i and j 

(including those at the junctions and their neighboring vertices), ^[ij] m ^2^ ^ s 
the sum over nearest neighbor bonds i and j on the skeletons between junctions, 
m ^2 * s the sum over bonds (ij) connecting the central point of the 
junction i and the nearest neighbor vertex j. The symbol tj in Sf ' is the unit 

(2) 

tangential vector of the bond i. The symbol in S 2 is the unit normal vector 
of the triangle i. 

S^ is defined only at the junctions. Note also that the bond [ij] in S^ 
includes the bonds connected to the central point of the junction. Si includes 
bonds that connect the vertices nearest neighbor to the junctions. The reason of 
this is for the sake of the in-plane elasticity at the junctions. When Si excludes 
such bonds, skeletons at the junctions freely move into the in-plane directions. 
The potential Si describes a one-dimensional interaction between vertices in the 
skeletons and a two-dimensional interaction only at the junctions. 

The partition function of the model is defined by 

,/ N 

Z= HdX lC xp(-S), S = Si+b 1 s£ ) +b 2 S?\ (2) 
J i=l 

where bi, 62 are the one-dimensional bending rigidity and the two-dimensional 
bending rigidity. In this paper, 62 is fixed to 62 = 5, and b\ is varied in the MC 
simulations. The center of the surface is fixed in the MC simulations, and this 
is denoted by / in f' Y[f=i dXj,. 

The canonical Metropolice MC technique is used to obtain the mean value of 
physical quantities. The vertices X of the skeleton are sequentially shifted so that 
X' = X+SX , where SX is randomly chosen in a small sphere. The new position X' 
is accepted with the probability Min[l, exp(— AS)], where AS = S (new)— S (old). 
The vertices of junctions (hexagons and pentagons) are shifted in two steps. 
Firstly 7 (or 6) vertices at one junction are shifted simultaneously at random, and 
the second step is a random translation and a random rotation of the hexagon 
(or the pentagon) as a rigid object. Those shifts are controlled by small numbers 
given at the beginning of the simulations to maintain about 50% acceptance 
rates. 

Four different sizes are assumed for surfaces in the MC simulations. The 
surface size (N,N s ,Nj,L) are (19362,10890,1212,4), (10242,5760,642,4), 
(5762,3240,362,4), and (2562,1440,162,4). In the simulations, two parameters 
are fixed: the total number of vertices L in the skeleton between junctions and 
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the two-dimensional bending rigidity 62 are fixed to 

L = 5, 6 2 = 5. (3) 

3 Results 

(a) 6i = 2.9, crumped (b) &i = 2.93, smooth 
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(c) the surface section of (a) (d) the surface section of (b) 



Fig. 2. Snapshots of surfaces of size (N,Ns,Nj,L) = 
(19362,10890, 1212, 4)obtained at (a) bi = 2.9 (crumpled) and (b) 61 = 2.93 
(smooth), (c) The surface section of (a), and (d) the surface section of (b). 



First, we show snapshots of surface obtained at b\ = 2.9 (crumpled phase) and 
6i = 2.93 (smooth phase) in FigsUfa) and|2Jb), respectively. The size of surface 
is (N,N s ,Nj,L) = (19362,10890,1212,4). The surface sections are shown in 
Figs|2Ic)|2Id). We clearly hnd that the surface is crumpled (smooth) at 61 =2.9 
(6i = 2.93). 

Figure Ofa) is plots of S%/N against 61 obtained on the surfaces of size 
(19362,10890,1212,4), (5762,3240,362,4), and (2562,1440,162,4). Because of 
the scale invariant property of the partition function Z in Eq.Q, S\/N is pre- 
dicted to be Sx/N = 3N/2(N - 1) ~ 1.5. We hnd from the hgure that the pre- 
dicted relation for Si/N is satisfied. The line connecting the data was obtained 
by multi-histogram reweighting technique. 
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(a) (b) 



Fig. 3. (a) Si/N against b%, and (b) X 2 against b\. b-i is fixed to 62 = 5. The 
curves are drawn by multi-histogram reweighting technique. 



Figure ffib) is the mean square size X 2 defined by 

i i 

where X is the center of the surface. The size X 2 of surfaces increases with 

increasing b\ and varies rapidly at intermediate b\ on the surface of 

(19362, 10890, 1212,4). This indicates the existence of the crumpling transition. 

The bending energy /Ny against bi is plotted in Fig. Ufa), where Ny = 
N—Nj is the total number of vetices where is defined. We find again a rapid 
varying of S% /Ny against b\ as the surface size increases. 

Figure 0Tb) shows the specific heat C„<i) defined by 

C s ^=^-((s^-(S^))\ (5) 

The curves drawn by multi-histogram reweighting technique indicate an anoma- 
lous behavior of C„(i). This can be considered as a sign of the phase transition. 
We expect that the transition is of first-order because of the sharp peak in C g (i> 

at the surface of (N, N s , Nj, L) = (19362, 10890, 1212, 4). The junctions play a 
nontrivial role in the phase transition, because we know that one-dimensional 
linear skeleton such as an elastic circle has no phase transition. The finite-size 
scaling analysis can be done with more extensive MC simulations including larger 
surfaces rather than those in this paper, and the first-order transition will be 
confirmed. 
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(a) ' th (b) ' b, 



Fig. 4. (a) S~ /Ny against bi, and (b) C <i> against bi. N v =N — Nj is the 

total number of vetices where is defined. 6 2 is fixed to b 2 = 5. The curves 
are drawn by multi-histogram reweighting technique. 



4 Summary and Conclusion 



We have investigated the phase structure of a skeleton model for membranes 
by using Monte Carlo simulation technique. Linear skeletons are joined with 
junctions and form a surface, which covers a sphere. Hamiltonian is a linear 
combination of the Gaussian bond potential, one-dimensional bending energy 
and the two-dimensional bending energy which is defined only at the 
junctions. Four different sized surfaces were used in the simulation. The length 
L of skeleton between junctions was fixed to L = 5 in those surfaces, and the 
two-dimensional bending rigidity b 2 was also fixed to 62 = 5 . 

We found that there are two distinct phases in the model; the smooth phase 
and the crumpled phase, which are characterized by small and large 
respectively, and also by small X 2 and large X 2 . Moreover, these two phases 
are separated by a phase transition, which was predicted to be of first-order on 
account of the anomalous behavior of the specific heat for 

Although we have obtained a sign of first-order transition in the skeleton 
model, more extensive MC simulations are necessary to confirm the result. The 
junctions are expected to play a non-trivial role in the phase transition, because 
we know that one-dimensional linear skeleton such as an elastic circle has no 
phase transition. Then, it is interesting to study the dependence of junction 
elasticity on the phase transition. A rigid junction model is also interesting. 
Rigid junction has an infinite in-plane and bending resistance against force, 
while the junction in this paper has finite elastic resistance. Phase structure of 
skeleton models can be clarified by further numerical studies. 
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